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Abstract 

In this paper the study of asymptotic stability of standing waves for a model of Schrodinger equation with 
spatially concentrated nonlinearity in dimension three, begun in [5], is continued. The nonlinearity studied is 
a power nonlinearity concentrated at the point a; = 0 obtained considering a contact (or S) interaction with 
strength a, which consists of a singular perturbation of the Laplacian described by a selfadjoint operator H^, 
and letting the strength a depend on the wavefunction in a prescribed way: iix = HaU, a = a{u). For power 
nonlinearities in the range (■;^) 1) there exist orbitally stable standing waves and the linearization around 
them admits two imaginary eigenvalues (absent in the range (0, -^) previously treated) which in principle 
could correspond to non decaying states, so preventing asymptotic relaxation towards an equilibrium orbit. 
This situation is usually treated requiring the validity of a nonlinear Fermi golden rule, which assures the 
presence of a dissipative term in the modulation equation ruling the complex amplitudes z{t), z{t) associated 
to the discrete part of the linearized spectrum. Here without the use of FGR it is proven that, in the 
range for a certain a* G (^, the dynamics near the orbit of a standing wave asymptotically 

relaxes in the following sense: consider an initial datum u(0) near the standing wave written in the form 
u{0) = uq = + e““+'’'‘’[(zo + + i{zo — Fq)'^ 2 \ + fo with zq small and /o small in energy and 

in a certain weighted space then the solution u{t) can be asymptotically decomposed as 

uit) = Too, as t ^ +(30, 

where Woo, koo > 0, G K, and tpoo and Too G L^(M^) , U{t) is the free Schrodinger group and 

IkoolU^ = as t ^ +00 . 

We stress the fact that in the present case and contrarily to the main results in the field, the admitted 
nonlinearity is L^-subcritical. 
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I. INTRODUCTION 


We continue here the analysis of a model of nonlinear Schrodinger equation with a concentrated 
nonlinearity in dimension three begun in [3]. We recall that such a model is defined by the equation 

= Ha(u)U. ( 1 ) 

where the nonlinear operator is a point interaction (or “delta potential”) in dimension three 

with a strength a = a{u) depending on the wavefunction u in a prescribed way. Here the nonlinearity 
is of power type and focusing. More precisely, the domain of is given by 

D{Hc,(u)) = G ; u{x) = ct){x) + qGo{x) with ^ G (2) 

^>o.| 

and the action is given by 

Hau{x) = —A(l){x) , X G \ {0} . (3) 

The nonlinearity, displayed in the boundary condition, means that the value at zero of the “regular 
part” (p of the element domain is related in a nonlinear way to the so called “charge” q of the same 
element domain, which is the coefhcent of the “singular part” Gq. According to Q, our choice for 
the function a{u) is 

a{u) = —u\q\‘^'^, z/, (T>0 

When fj = 0 one obtains the well known contact interaction (of strength a = —i'), which, due to 
the sign, is a so called attractive 5 interaction (see ID)- When (T 7 ^ 0 the nonlinearity is in some 
sense acting at the single point zero, coinciding with the location of the singularity of the contact 
interaction. This is the origin of the denomination of concentrated nonlinearity (see [DEI IE]). 
The above model can be derived from a standard nonlinear Schrodinger equationquation with a 
inhomogeneous nonlinearity, i.e. with a nonlinearity space dependent, and shrinking at a point in a 
suitable scaling limit. This derivation is rigorously treated in |in) for the one dimensional case and 
in the forthcoming paper m for the present three dimensional case. The problem 0 describes a 
Hamiltonian system for which global well posedness holds in the nonlinearity range a G (0,1). More 
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precisely we endow the space L^(M^,C) ~ (assuming the usual identification 

z ~ with the symplectic form 



uv dx 


[ (JRuQ'n — '^v?R.u)dx 
Jr3 


{u2Vl — UlV2)dx 


(4) 


Of course M)©ff^(M^, M) is a symplectic submanifold, and we associate to (Q the Hamiltonian 

functional coinciding with the total conserved energy associated to the evolution equation ([^, that 
is given by 


E{u{t)) 


y^ = 4> + qGoev. 


(5) 


where V is the domain of finite energy states 


V = {ue L2(m3) . + qGo{x), with 0 € V(j) G L‘^{R^), q E C}, (6) 


which is a Hilbert space endowed with the norm 


klly — II^</'IIl2 + \q\" 


(7) 


Note that for a generic element u of the form domain the charge q and its regular part 
independent of each other. 

Correspondingly, the NLSE 0 can be rephrased in the hamiltonian form 

du ^ , 

- = J E (u(t)) . 


are 


( 8 ) 


0 1 \ 

where J = \ is the standard symplectic matrix. In [3] the case of power a E (0, was 

,-10/ 

studied, showing existence of nonlinear bound states, orbital stability and relaxation to (relative) 
equilibrium asymptotically in time. For such values of the power nonlinearity a the asymptotic 
analysis is simplified by the fact that linearization around standing waves has no eigenvalues apart 
from the zero eigenvalue, always existing due to the gauge 17(1) symmetry. For a E (:/^)l) the 
linearization around a standing wave admits pure imaginary eigenvalues which correspond to 
existence of neutral oscillation in the linearized dynamics. If these neutral oscillations persist as 
invariant tori in the phase space of the complete nonlinear system, relaxation to an asymptotic 
equilibrium standing wave is precluded. The analysis of so called asymptotic stability of solitons 
started at the beginning of the nineties with the studies of Soffer and Weinstein (MM) on the 
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NLS equation with an external potential in dimension three and of Buslaev and Perelman on the 
translation invariant NLS in dimension one ([318]). 

In our model the analysis goes as follows. We consider the nonlinear evolution problem 



u{0) =Uo £ D{Ha), 


(9) 


In [3j it is shown the existence of a solitary wave manifold for Q 


M = { = 


4:TTU 


47r|x| 


: w > 0 


( 10 ) 


and it is shown that these solitary waves, which belong to D{Ha), are orbitally stable in the same 
range (a G (0,1)) where global well posedness of equation Q is guaranteed. 

One aims at proving that an initial datum near this solitary manifold relaxes for t —)> +oo to the 
solitary manifold itself. 

The solitary manifold turns out to be a symplectic two dimensional submanifold. 

Now writing u = e*‘^*(<I>tj + R), we obtain that R satisfies the first order the linearized canonical 
system 


dR 

dt 


J 

Hai ® LO 0 

R = 

0 

L 2 


0 Ha2 + w 


_-Li 

0 


R 


= LR , 


( 11 ) 


and Lj = Ha^ + uj for j = 1, 2, ai = -(2cr + 1)^ and a 2 = 

As recalled before (see |3|, Section IV A, and Appendix VI A of this paper) in the case a £ {l/y/2, 1) 
the discrete spectrum of L consists of the eigenvalue 0 with algebraic multiplicity 2 and of two purely 
imaginary eigenvalues with 


^ = 2cr\/1 — 


( 12 ) 


and corresponding eigenvectors 'L and T* (see Appendix VI C of this paper). As a consequence, the 
domain of the operator L can be decomposed in three symplectic subspaces, more precisely 

D{L) = 

where X^, X^, and are respectively the generalized kernel of L, the eigenspace corresponding to 
the L eigenfunctions T and 'k*, and the subspace associated to the absolutely continues spectrum 
of L. In particular X^ is generated by the tangent vectors to the symplectic solitary submanifold. 
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The corresponding symplectic projection operators from onto X^, and X'^ are denoted by 

P®, P^, respectively (see appendix C for explicit representation .) A further step in the analysis 
consists in decomposing the solution of equation ([^ in the sum of a soliton-like part 
and a fluctuating part x(t, x), introducing the Ansatz 


with 


u{t,x) = + x{t,x)) , 


0(t) = / Lo{s)ds + ■y{t), 

Jo 


and 


xit,x) = z{t)^{t,x) + + f{t,x) = 'tjj{t,x) + / G 


(13) 


(14) 


(15) 


with ui{t), 7 (t), z{t) and f{t,x) up to now unprejudiced. Notice that the parameters are now time 
dependent, as well as the functions T and T: this is due to the fact that 'h depends on t through L 
and then through w. 

The goal now is to show that the fluctuating part is decaying (4'-4;'* component) or dispersing (/- 
component), to provide convergence of parameters to possibly unknown asymptotic values; 

all of this finally gives relaxation to the solitary manifold. To this end one has to fix the above under¬ 
determined representation, and then obtain equations for 7 (t) and x(f,x). This is achieved by 
requiring the fluctuation x to be symplectically orthogonal to the solitary manifold At (or, equiva¬ 
lently, orthogonal to the generalized kernel Ng{L)) for every time t > 0. This procedure yields the so 
called “modulation equations” for w(f), 7 (t) and x(t,x). Moreover, exploiting further orthogonality 
relations between one also gets equations for the coefficient z and the dispersive term / 

in X (See Theorem II. 1 in Section |l^ for more details). At this point two problems occur. The first is 
that the modulation equation for the fluctuating component x, due to the introduction of the time 
dependent Ansatz, contains a non autonomous linear part. So, the use of dispersive estimates to 
show the wanted decay of this component requires to preliminarily “freeze” the dynamics at a certain 
fixed time T; the subsequent step is to show the uniform character of the estimates in T. The frozen 
equations are written in Section |IIj In the same section the leading terms in the modulation equa¬ 
tions are identified, and the estimates on the remainder terms are displayed. This is however not yet 
sufficient to get rid of the complex oscillation described by ^ and z*. So the modulation equations 
are rewritten in a canonical way by means of a Poincare normal form pushed to the third order 
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which preserves estimates on the remainders. In particular, the transformed oscillating component 
zi = zi(z,z) satisfies the equation 


zi = i^zi + iK\zi\^zi + Zr 


(16) 


The asymptotic behavior of the solution of the zi component depends on the coefficient iK and more 
precisely on the sign of its real part. 

This is the point where nonlinear Fermi golden rule (FGR) enters the game. It turns out in relevant 
examples that if a resonance condition between a higher harmonic (a multiple) of and the continu¬ 
ous spectrum of linearization is satisfied (the FGR), then ‘R{iK) is strictly negative; this gives decay 
of the oscillating modes of the linearization (see the seminal paper cited above and moreover |24] , |25) 
and m)- The decay is ultimately due to coupling of oscillating modes with the continuous spectrum 
assured by FGR, and consequent drift of energy from the discrete component to the continuous one; 
so the mechanism is dissipation by dispersion. Furthermore, exploiting the hamiltonian structure of 
the system it can be shown that the above situation is in some sense generic ([I3l[13j,l5] and reference 
therein). 

In the present model things go as follows. In the first place the second harmonic 2i^ of the 
discrete eigenvalue lies in the interior of the continuous spectrum if the nonlinearity power a G 
^(73’ ^ 27 ^^’ then denote by N 2 {q, q) the quadratic terms coming from the Taylor expansion of 

_ { Ti(u;o) \ 

the nonlinearity (see formula (20)), and by 'I'(a;o) = I I the eigenfunction of the linearized 

y 'I'2(wo) j 

operator associated to (given explicitely in Appendix VI C, Proposition VI 4). The nonlinear 
Fermi golden rule (FGR) adapted to our case is the following non-degeneracy condition: 




(17) 


where 4'_|_(2z^o) is the generalized eigenfunction associated to -|-2z^o (for the explicit representation 
see Appendix D). 


This nondegeneracy condition should imply ^(iK) < 0 in (16). 

As a matter of fact, thanks to the explicit character of our model, we are able to directly verify, 
without use of the nondegeneracy condition, that ^{iK) < 0 (that gives dissipation in ([I^) holds 
for any a in the range for ^ certain a* G Section HID). Moreover, we 

have numerical evidence that this is true on the whole interval With these premises, 

we eventually prove the following result. 
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Theorem (Asymptotic stability in the case of purely imaginary eigenvalues) Assume that 
u £ (7(1^+, V) is a solution to equation Q with a power nonlinearity (see Q) given by u G 
for a certain a* G Moreover, suppose that the initial datum is close to a standing wave 

of in the sense that 

u(0) =uo = + i{zo - ^)^2] + fo ^ V n Li{R^), 

with loq > 0, 7 o, zq G M, and fo G L^(M^) n -L^(M^) and 

|zo| < and ||/o||li, < ce^^^ 

where c, e > 0. 

Then, provided e is sufficiently small, the solution u(t) can be asymptotically decomposed as 
u(t) = as t -£ +oo, 

where ujoo, e^oo > 0, 6i G M, and V’oo ?’oo £ T^(M^) such that 

lkoo||L 2 = 0(t“^/^) as t —> + 00 , 


in L2 (m3). 

Some remarks are in order. 

We stress again that the range of the admitted nonlinearities a implies that is in the essential 
spectrum of the linearized operator. Notice that ~ 0.96, and that when cr —>■ 1 the discrete 

eigenvalues of linearization, ±i^, collide at zero. So, in order to cover a larger range of values of a 
one has to consider harmonics higher than the second and consider normal form pushed to order 
higher than the third. As shown in [3] the standing waves in the case cr = 1 are unstable by blow-up. 
We recall that u = is a threshold resonance for the linearization. Its presence does not allow 
enough decay of the linearized dynamics to obtain relaxation to solitary manifold. 

The asymptotic stability result is achieved following the outline of [9] and |18j and using the machin¬ 
ery already developed in [3]. In particular, in [18] the same problem for the analogous one-dimensional 
model is studied. 

Nevertheless, the three-dimensional case presents some differences. The first one is that the concen¬ 
trated nonlinearity imposes to develop the analysis at the form level. This means that the estimates 
on the evolution of the initial data are more delicate. 
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The second main difference is the faster decay of the propagator of the free Laplacian. This allows to 
develop the the analysis using just the structural weight rc = 1 + |^ which arises from the dispersive 
estimate instead of introducing new weighted spaces as done in the one-dimensional case treated in 

mm- 

Finally, the eigenfunctions associated to the purely imaginary eigenvalues do not have oscillating 
terms occurring in the one-dimensional case but they exponentially decrease as |x| —)> -|-oo. This fact 
will be useful in order to get the decay in time of the radiation term. 

A last comment of general nature is in order. As in the one dimensional case studied by Buslaev, 
Komech, Kopylova, and Stuart in [6] and Komech, Kopylova, and Stuart in [18], and the three 
dimensional model analyzed in |3], the analysis of a specific model allows to obtain asymptotic 
stability of standing waves without a priori assumptions. In particular the nonlinearity is fixed, of 
power type and subcritical, in the sense that it falls in the range of global well posedness of the 
equation (see m for a different example where asymptotic stability is proven in subcritical regime); 
no spectral assumptions are needed; no smallness of initial data is required, in the sense that we give 
results for every standing wave of the model and initial data near the family of standing waves. 

As a final remark, while Komech, Kopylova, and Stuart in [T5] find a link between the Fermi Golden 
Rule and the decay of normal modes of linearization, here such decay is directly verified. This fact 
seems to indicate that some of these assumptions or hypotheses are in fact unnecessary when enough 
information about the model is known. 

For the sake of completeness, in the course of the paper we will repeat proofs requiring some modifi¬ 
cations because of the facts mentioned above; on the contrary, where the arguments hold unchanged, 
just a reference will be given. Moreover in the appendices we give information about linearization 
operator recalling useful material from [3| and giving detailed properties of discrete and continuous 
eigenfunctions. 

Acknowledgments The authors are grateful to Scipio Cuccagna, Gianfausto Dell’Antonio, Alexan¬ 
der Komech and Galina Perelman for several discussions and correspondence. R.A. and D.N. are 
partially supported by the FIRB 2012 grant Dispersive Dynamics: Fourier Analysis and Variational 
Methods, reference number RBFR12MXPO. R.A. is member of the INdAM Gruppo Nazionale per 
I’Analisi Matematica, la Probabilita e le loro Applicazioni; D. N. is member of the INdAM Gruppo 
Nazionale per la Fisica Matematica. 


9 


II. MODULATION EQUATIONS 


We recall that the operators we are dealing with have different domains while the forms associated 
share the common domain V, introduced in Q. It proves convenient to describe V in an alternative 
way: fixed an arbitrary A > 0 and denoted G\{x) := one finds 

V = {u = (f)\ + qG\, with(/)A G q € C] . 


The arbitrary positive parameter A is customary in the description of a point interaction and more¬ 
over allows the regular part of the element’s domain to be in the Sobolev space instead of the 
corresponding homogeneous counterpart, which is often useful. Everything is in fact independent on 
the choice of A . 

Correspondingly, we will perform computations at the form level. In order to do that let us recall 
that the variational formulation of equation 0 is 

= Qa{u{t),v) Vu E U. (18) 


where, given u,v £V with u = 4>u,x + QuGx, v = (f)y^x + QvGx, 

Qa{u,v) = iV(j)u,\,^4>v,x)L2 + quQv (a + j - Xqu{Gx, (t>v,x)- Xqv{4>u,x,Gx)L^ 


Note that the equation (18) makes sense because V is independent of the positive parameter A and 


it is a Hilbert space with the norm 


u 


|y = l|V(/)A|li2 + kl^ yuGV. 


In order to inspect the asymptotic stability of equation 0 it is useful to represent the solution u by 


the aid of the Ansatz (13) together with deflations ( |14| ) and (15). From now on we always refer to 
such formulas. 

Hence, we want to construct a solution of equation 0 close at each time to a solitary wave. Let 
us notice that the solitary wave does not need to be the same at every time, which means that the 
parameters uj{t) and 0(f) are free to vary in time. 

Exactly as in the case a E ^0, (see [3], Section V) the function y solves 


■ dx 

i— 

dt 


it),v] = Qa,Linixit),v) +7(f)($a;(t) +X{t),v)L 2 + 

J L2 


(19) 
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+Uj{t) ( -i 




U)(t) 


doj 


,v) + N{q^{t),q^), 

/ L2 


for all V \ 

Here Qa,Lin is the quadratic form of the linearization operator acting as 


Qa,Lin{x, v) = (pv )“ -^^{Qx^v) “ CT^?flq^?flqy + u;(x, v)l2 , 

and the nonlinear remainder N[q^, q^) is given by 

N{(lxi + 2 ^( 20 - + l)\quj\^'''Rqx^qv + u\q^\‘^‘^^q^^qy + i^\qui\‘^''?fl{quj^), 

( 20 ) 

J_ 

where (see section II B in |3]), ■ 

Since u{t), 7 (t), and unknown and the propagator grows in time along the directions of 

the generalized kernel of the operator L, the idea is to get a determined system requiring the function 
xit) to be orthogonal to the generalized kernel of L at any time t > 0. Hence, one obtains that w, 
7 , z, and / must solve the following system of equations. 


Theorem II.1. (Modulation equations) If xit) is a solution of equation (19) such that P^xit) = 0 
for all t > 0 and u){t) and x{t) are continuously differentiable in time, then u and 7 are solutions of 


u = 


+ x) 




3? {JN 


7 = 


^-xV 




( 21 ) 


( 22 ) 


Furthermore, z and f satisfy 


{'Ii,J^)l 2 {z -ifz) = iR.{J N {q^)qjq,) + uj 


uf) -(^.Jf 

d^Jh^ \duj / l2. 


+ 'iix,^)p, (23) 




+ {87TVXP^JN{q^)Gx,qvGx)L^, 


( 24 ) 


for all V ^ V. 
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Proof. Equations (21) and (22) can be proved with the same argument exploited in the case a G 


0, (see |3], Theorems V.3 and V.6). 


Equation (23) can be obtained taking v = JT as test function and noting that 


= 0 , 


Einally, equation (|24|) follows taking the projection onto the continuous spectrum of both sides 

□ 


of equation (19) and recalling that / G 


A. Frozen spectral decomposition 

The goal of this subsection is to get an autonomous linearized equation for the component /. 
According to (0, having fixed T > 0 we define 


Lx = L(^ujt) '■= ~J 


^ai{T) + 0 

0 Ha2{T)+>^T 


where Q!i(T) = —{2a + 1)^^^^, a 2 = — and ujt = w(T). Then for any t G [0, T] one can 

decompose f{t) G X'^ = X^{t) as 

f = g + h with g G X!}^ = X^(B X^, hG 

where the subscript T means that the time is fixed at t = T. We recall that and P^ are the 
symplectic projections on the kernel and nonzero eigenvalues of the linearization (see Appendix C). 
Correspondingly we denote by P^ and Pf the projections on the analogous subspaces frozen at time 
T. Finally denote P^ = P^ + P^. 

Then for the quadratic form 


Ql{u,v) - Qlj,(u,u) = ^ '''^ -?R.{Tquqv) - (wr - ^){Ju,v)l2, 

47r 
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for all tt, V G V, where 


T = 


0 -1 
2a+ 1 0 


Hence, observing that = 0, the equation (24) for / is equivalent to 


df 


dt' 


1,2 




-77,v] =QLT{f,v) + {{uj-UT)Jf + d}——'il} + 'yP‘'Jx,v] + 


duj 


( 25 ) 


L2 


+ (sttVX ^^ Tqf + P^^JNiq^)^ Gx, q,Gx^ 


L2 


for all V € V. 

Since our dispersive estimate holds only on the continuous spectral subspace, we need to prove that 
it is enough to estimate the symplectic projection of x(i) onto that subspace. This is stated in the 
following lemma where we denote by TZ{a) a bounded continuous real valued functions vanishing as 
a —)■ 0, and 


Pl(w) — P{\\U] - Wo||c'0([o,T]))- 

Along the paper, with a slight abuse, we shall use the symbol P(a, b) to denote bounded continuous 
real valued functions vanishing as a, b ^ 0. 

Lemma II.2. If \oj — ujt\ is small enough, then the function g can be estimated in terms of h as 
follows: 




The last lemma can be proved following the proof of Lemma 3.2 in |18] . 

As a consequence, one can apply the operator Pf, to both sides of the equation for / 


^ ~ + { Pf 


dP^ 

{u - 0 Jt)J f + dj— —V’ + fP^JX 
duj 


L2 


and obtain 
+ (26) 


+ (^TrVXPf^ ^f^ Tqg + P^JN{q^)^ Ga, q.Gx 


L2 


for any v &V. 
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B. Asymptotic expansion of dynamics 


In order to prove the asymptotic stability of the ground state we need to show that for large 
times 2 ; and h are small. For this purpose, in this section we expand the inhomogeneous terms in the 


modulation equations, then rewrite the equations of cj, 7 , 2 : and h as in (29), (30), (31) and (36), 
and for each equation we estimate the error terms. 

In what follows we denote 


= qiPi + <12P2, Vp, g G 




1 /( 2 < 7 ) 


With an abuse of notation we denote by J charge of the function 




As a preliminary step, we expand the nonlinear part of the equation (19) N{q-^ as 


N{q^) = N2{q^) + N^{q^)+NR{q^), 


(27) 


where N 2 and are the quadratic and cubic terms in q^ respectively, while Nr is the remainder. 
Exploiting the Taylor expansion of the function F{t) = around \quj\‘^, one gets 

^{N2{qx)^) = ^{{cr\qu;\'^^''~^^\qx\‘^qLo + 2a\qu;\‘^^'^~^\q^,q^)q^ + 2{a - l)cr\q^\‘^^‘^~‘^\q^,q^fqu;)^), 


and 


^{Fl 3 {q^)q^) = ^^\qy^\‘^q^ + 2{a - l)a\q^\‘^^‘^ “^^q^, Qxf Qx+ 


+2{a - i)cr|g^p('^ ‘^\qu,,qx)\qx?‘<i^ + - 2)((t - i)cr|g^p('^ '^\qu:,qxfq^)qv)i 

for any q^ G C. For later convenience, let us define the following symmetric forms 

N2iqi,q2) = cr\quj\‘^^‘'~^Hqi, q2)qu} + cr\qoj\‘^’''"~^'>[iqoj, qi)q2 + {qu},q2)qi] + 


+2{a - i)cr|g^|^('^ ‘^\qu,qi){q^,q2)qu. 


and 


1 3 ^ 3 

^^3(91,92,93) = ^ (9i,9i)9fc + 3(0--l)o-|9a;|^^‘^ (9^,9i)(9c^,9i)9fc+ 
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1 ^ 4 

+ -{a - {qoj,qi){qj,qk)quj + -{a -2){a - l)(T\qu:\‘^^'"~^\qu,, qi){quj, q 2 ){quj, q^)quj- 

2,ji,fc=l 

In order to prove the asymptotic stability result, we shall prove in Section 2.4, the following asymp¬ 
totics 

ll/Wlli”_i ~ z{t)r^t~^, ||V’(t)||i/~ ( 28 ) 

w ^ 

as t —)■ 4 - 00 . 


Remark II.3. As in m, the first step in proving these expected asymptotics is to separate leading 


terms and remainders in the right hand sides of the modulation equations (21) - (23), (26). Basically, 
in the next subsections, we will expand the expression for oj, 7 , and z up to and including the terms 
of order t~‘^l’^, and for h up to and including t~^. 


Remark II.4. Note that since the nonlinearity depends only on the charges the same holds for its 
Taylor expansion. 


1. Equation for uj 


Substituting the expansion for the nonlinear part N given in (27) in equation (21) and considering 


the asymptotics (28) one gets 


1 


w = —ifi{{JN2{q^)+2JN2{q^jj,qf) -4 JN^{q^j,))qu) + { V', 


A 2 


d^ 

doj 


L 2 


WN2{q^)qu>) + ^R, 


where A = and the remainder VLr is estimated by 


\^r\ < Tliuj, |^;| + ||/||l°°_ 

lit ^ til 


Recalling that = z'k -|- z'k*, one can rewrite the previous equation for to as 
UJ = ^l2oz‘^ + ^iizz ^02z'^ + -|- id,2iz'^z ^i2zz‘^ -|- 4“ ^{qfi ^ 10 ) ^iqfi ^ 01 ) T Ilij, (29) 


where the rijj’s and the are suitable coefficients, while fl/j is a remainder term. 
Remark II.5. Since the second component of the vector equals 0, one has 

fill = 2^Re{JN2{qi;,)q^) = 0 . 

This fact will turn out to be useful in writing the canonical form of the modulation equations. 
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2. Equation for 7 


As in the previous subsection the equation for 7 (22) can expanded as 


1 1 / 

7 = —^{{JN2{q^) + 2JN2{q^,qf) + JN^(q,b))qjd^^) + ^ 


where the remainder is estimated by 

\rn\<n{u;M + \\f\\L- ,){\z\^ + \\f\\L- ,f- 

As before, the equation for 7 shall be written in the form 


7 — ^2oz'^ + ^iizzTq2z'^ + Fso^^ + T2iz‘^z + Ti2Zz‘^ + Fqs^;^ + z{qf, F^^q) + z{qf, Fgj^) + F/j. ( 30 ) 


Remark II.6. In this case Fn does not vanish as in equation ( 29 ). 


3. Equation for z 


Exploiting the results of the previous subsections, equation ( 23 ) can be expanded as 


z-i^z = \'R{JN2{q^)qf) + ]:'R{{JN2{q^) + J N'i{q^))qj^) + 


where k = —('F, J'F)^2 and 


ID ID 


Mimicking the notation employed in ( 29 ) and ( 30 ), the previous equation can be written in the form 


i = if,z + Z20Z + Z11ZZ+Z02Z + Z^qz'^ + Z21Z z + Z12ZZ + Zosz'^ + zft{qfZ[Q) + zft{qfZQ-^^) + Zr, ( 31 ) 


and it turns out that 


Zii = ^5R(JA'2(9f, ^f*)?’^)) Z 20 = \fl{JN2{q^)q^i)., Zq 2 = \^{JN2{q^*)q'^)-, 

■^21 = ^'^{JN^{q^*,q^,q^)^) + ^ J^)l2 ^{J^2{q^)qj<t>J)+ 

-{^\'I>)L2^{JN2{q^)q^)-2\ml,^{JN2{q<i.^,q<i.)q^ 

j. . ^ j, . 


Z(o — 2 


JN2{q^*m) 

K ■ 


did 

r/l _ 

^01 — ^ K 


( 32 ) 
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4- Equation for h 


In order to expand asymptotically the equation (26) for h, the following remark will be useful 


Remark II.7. For any f G the following holds 


P^P^f = P^{I - p^)f = P^{P^ + P^- P'^)f = P4pf + P^{P!}, - P'^)f 


Let us denote 


p{t) =uj{t) -wr + 7(t), 


(33) 


then equation (26) can be rewritten as 


—, ^)i/2 + {8TrP^JN2{qp)G\, qvGx)]^2P 


+ P 


uj- 


dP^ 

duj 


^ + pJg + 7(R^ - P^)Jf 


,V I + 

p 


+ (^SttVXP^ + P"JN{q^) - JN^iqp)^ Gx,q,Gx^ ^ , 

for any v € V. 

Denote 


H'r = Pt 


(jpc 

+ iP^Ji’ + pJg + UPt - P‘')Jf 


and 


H'^ = SttVXP^ + P^JN{q^) - JN2{qp)^ Gx- 

We recall that by 11^ we denote (see Appendix C) the projections onto the branches C± of the 
continuous spectrum separately. Analogously we denote by 11^ the corresponding projections of the 
linear generator Lt frozen at time T. The following lemma is useful in the following. 


Lemma II.8. There exists a constant G > 0 such that for each h G holds 


\\[p-j-i{U+-U:f)]h\\^, <G\\h\\L^_^. 

" " -^w W 
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The proof is in Appendix VID for any t > 0. Finally, let us define 


LM{t) = Lt + - n^), 


( 34 ) 


then the previous equation becomes 


dh 


—, u ) — qvG\)i2 + {Hu, v)i2 + {H^, q^G\)i2, (35) 


L2 


for any v € V, where we denoted 


Hr = H'j^ + p[P^J - i(n+ - ny)]/i. 


Finally, let us expand the second summand in the right hand side of (35), and get 


dh 


— ,uj — QL^j{h,v) + {z H 20 + zzHii + z Ho 2 )qv + {Hr,v)r 2 + {HR,qyG\)i 2 , (36) 


for any v € V, where 


H 20 = {87rVXP^JN2{q<,)Gx,Gx)L^, 

Hii = 2{8'Kxf\P^JN2{qm,qm*)G\,Gx)R2, 

Ho 2 = {8nVXP^JN2{q^*)Gx,Gx)L^. 


Thanks to the estimates done for the other equations and Lemma II.8 one can estimate the remain¬ 
ders in the following way: 

WP'rWlI < G (|2:|(|cj| -h \^\)+ni{uj){\uj-UT\ + |7|||/||l-_i)) < 

< Tli{uj, |z| -t- ||/||ioo_ ) (\zf -I- |2;|||/||£,oo_ -I- ||/|||oo + Iw - wr|||/||L°°_i) , 

ii; ^ \ w w~^ / 


hence 


II^rIIli, < H\{^, kl + \\f\\L'^_-^) (1 




L“_i + 


J^OO ^ “1- \0J - U)'J' 


(37) 


and 


\Hr\\l]^ < 77i(w, |z| -I- ll/ll^oo ) (\zf p k|||/||L°°_^ + ||/||i°° + |w - cjtKI^P + II/I|l°°_i)) ■ (38) 

^ It; \ w w~^ w ^ / 
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Remark II.9. In the same way one directly expands the equation for the function f getting 

+ 2;^Fo2)(7« + (-Fr, + (F^, g^GA)L2, (39) 

for any v G V, where 

F 20 = {87Ty/\JN2iqis,)G\,G\)i2, 

Fn = 2{8TT^/\JN2{qm,q^!*)G\,G\)i2, 

Fo 2 = {8Tl^/\JN2{q^*)Gx,Gx)L2. 

and 

— (iP^ 

Fr = u—fj + + j{p^ - p^)jf, 

F'f, = 87^^/A { ^~^f^ Tqf + PnN{q^) - JN2{q^)^ Ga- 

Furthermore, the norms of the remainders Fr and F'f^ can be estimated by the corresponding 
norms of the remainders Hr and H'f^. 

III. CANONICAL FORM OF THE EQUATIONS 

In this section we would like to use the technique of normal coordinates in order to transform the 
modulation equations for oj, 7 , z, and /i to a simpler canonical form. We will also try to keep the 
estimates of the remainders as much close as possible to the original ones. 


A. Canonical form of the equation for h 


Our goal is to exploit a change of variable in such a way that the function h is mapped in a new 
function decaying in time at least as For this purpose one could expand h as 


h = hi + k + ki, 


(40) 


where 


k = 020-2^ + anzz + 002 ^^, 

with some coefficients aij = aij{x,uj) such that 07 = aiji, and 

ki = —exp LMis)ds^ /c( 0 ), 
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where the operator Lm was defined in (34). 


Note that hi(0) = /i(0), since A:i(0) = — A:(0). 


Proposition III.l. There exist aij G for i, j = 0, 1, 2, such that the equation for hi has 

the form 


,v j — Qlm {hi,v) + {Hr, v)i 2 + {Hr, qyG\)R 2 , 


(41) 


for all V ^ V, where Hr = Hr + Hr with 

Hr = - 


^ /da^^2 _|_ j _|_ {2a20Z + aiiz){z - if,Tz)+ 

\ dio doj doj ' 


(42) 


+(aii2; + 2a2oz){z + if,Tz) - /o(n^ - n^)A:] . 


Proof. The thesis is proved substituting (40) into (35) and equating the coefficients of the quadratic 
powers of 2 ; which leads to the system 


QLt{<^20,v) +^{H20^) — {Hf,Ta20,v)L2 = 0 
Qlj, (an, u) + 3f?(i7ii^) = 0 
Qlt{^02,v) + 3^(-ffo2^) + {Hf,Tao2,v)L^ = 0 

for all V € V. The former system admits the solution 

an = —LjfHii 

020 = —{Lt — 2 fCr — H)~^H2o 


002 = 002 = —{Lt + 2*Cr — 0) H^ 


02 


(43) 


□ 


Remark III.2. From the explicit structure of the remainder Hr it follows that it still satisfies 


estimate (37). 


We will need to apply the next lemma which can be proved as Proposition 2.3 in 


Lemma III.3. If cr € ^ ^ HL\,, then there exists some constant C > 0 such that 

for any t > 0 

\\e-^Tt^LT + 2iiT - 0)-'P#/||l-_, < ^(1 + t)-^/^\\f\\Li. 
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Remark III.4. Let us note that 


h = P^h = Ptfhi + P^k + P^ki, 

hence, in order to estimate the decay of ||/i||^oo_ ^ it suffices to estimate the decay of 

w 1 

\\Pfki\\L°°_,- 

W ^ W W 

B. Canonical form of the equation for oj 

Since = 0, we can exploit the method by Buslaev and Sulem in [9], Proposition 4.1 and get 
the following proposition. 

Proposition III.5. There exist coefficients bij = bij{io), with i, j = 0, 1, 2, 3, and vector functions 
bT = bT{x,u), with i, j = 0, 1, such that function 

UJI =UJ + b20Z^ + buzz + bo 2 z‘^ + bsoz^ + 621 -^^^ + b^zz'^ + 603 ^^+ 


+z{f, ^io)l2 + ^(/) ^oi)l2, 
solves a differential equation of the form 

= LIr, 

for some remainder LIr. 


Proof. Substituting the equations (29), (31), and (39) into the derivative with respect to time of the 
expression for uji and equating the coefficients of z^, zz, z‘^, z, and z one gets the following system 

+ 2^^620 = 0 

^02 — 2i^6o2 = 0 

^30 + 3^^630 + 2 ^ 20^20 + ^{F2oqbff) = 0 
^03 — 3i^6o3 + ‘2'Zo2bo2 + ^{Fo2qb^) = 0 
^21 + ifb 2 i + 2Z11620 + 2^20^02 + + -^20^^) = 0 

f2l2 — ifbi 2 + 2Z11602 + 2^20^20 + + -^ 209 b'jg) = 0 

(g/, O'lo) + if if, b[Q)L2 + QlU, 6'io) = 0 

(«/> ^ 01 ) + ^oi)l2 + Ql(/, 6'oi) = 0 


The last two equations of this system can be solved in a way similar to the ones system (43), and 
the proof follows. □ 
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Remark III.6. From the proof of the previous proposition it also follows that the remainder H/j can 
be estimated as Qr, namely 


< '^(w, | 2 ;| + \\f\\L^ ){\z\‘^ + 

lit ^ 'll! 


In the next lemma we prove a uniform bound for \ujt — uj\ on the interval [0, T], For later convenience 
let us denote 


TZoioo, \z\ + II fllroo ) = TZ I max \ojt — w|, max {\z\ + 

\0<t<T 0<t<T 


Remark III.7. Let us note that |u;| < |a;o| + Icjq — ^t\ + |a; — ojtI, then 



max Tl{u, 
0<t<T 



( max \ujt - a;|, max (|z| + ||/||l°o ^) 
\0<t<T 0<t<T 


The next lemma can be proved as in Section 3.5 of jl8j . 


Lemma III.8. For any t G [0, T] we have 


|wr - w| < 7^2(w, |z| + ||/||l°°_j) 


Jt 




+ (kT| + ||/T||L» J^ + (k| + ||/||L-,)^ 

lit lit 


C. Canonical form of the equation for 7 


Equations (30) for 7 and (29) for w differ just because in general Fn 7 ^ 0. But we can perform 


the same change of variable in the previous subsection, namely 

7 i = 7 + d2oz‘^ + do2z‘^ + dsoz^ + + di2zz‘^ + dosz^ + z{f, dio)L2 + z{f, 41)7,2, 

for some suitable coefficients dij = dij{uj), with i, j = 0, 1, 2, 3, and vector functions rfh = d[j{x,uj), 
with i, j = 0, 1. Then the function 71 solves the differential equation 

71 = rii(a;)zz + ri{, 

for some remainder F^j, which can be estimated as Tr, i.e. 


\fR\<n{u;,\z\ + \\f\\L-,){\z\^ + \\f\\L-J^ 

in J- in J- 
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D. Canonical form of the equation for z 


Exploiting the change of variable (40) used to obtain the canonical form of equation (35) for h, 


one can prove the following proposition. 

Proposition III.9. There exist coefficients Cij = Cij{uj), with i, j = 0, 1, 2, 3, such that function 

Z 1 = Z + C20Z^ + CnZZ + C02Z^ + C 30 Z^ + C 2 lZ^Z + C03Z^, 

solves a differential equation of the form 


zi = if,zi + iK\zi\‘^zi + Zr, 


(44) 


where 


* " " * -72 2 


iK — Z21 + .^21 + ^■^20-2^11 ~ ~ 


with the coefficient Zij, i, j = 0,1,3, defined in (32), and 

= (5 + Pf^i + Pr^i) + (5 + Pj^hi + Pfki, Zqi)z + Zr. 


The proof is a matter of calculation, but we give it explicitly to stress the role of the functions Uij, 
h j = 0, 1, 2. 


Proof. Substituting (40) in the equation (31) the differential equation for z becomes 


z — ifz + .^20-2^^ T Z\iZZ + Zq 2 Z^ + Z 3 qZ^ + Z 2 lZ^Z + Zy 2 ZZ^ + Zq 3 Z^-\- (^5) 


+-^30^^ “I" ^ “b T -^03^^ P ^R^ 


where 


Z'o = ^{qaA)^ 


Z '3 = ^{qaA)^ 


Z'^^=^{qaA)+^{<laA)^ 


Z[^ = AA+^{<laA)^ 
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and the remainder Zr is as in the statement of the proposition. 


Inserting equation (45) into the time derivative of the expression for zi and equating the coefficients 
of zz, z^, zz‘^, and z^ one obtains the system 


*^C20 + -^20 = 0 
—iicii + Zii = 0 

—3i^co2 + ^02 = 0 

2iCc30 + -^30 + -^30 + 2 c20^20 + C 11 Z 20 = 0 

Z12 + Z'^2 + 2C20-^20 + C.\\{Z\\ + Z02) + 2co2-Z'ii — 2z^Ci2 = 0 

—4i^C03 + ^03 + ^03 + Cll^02 = 0 


The theorem follows from the fact the the above system is solvable and in particular 




and 



□ 


Remark III.10. For later convenience let us note that, since Z 21 , Z 20 , Zu, and Zq 2 are purely 
imaginary, one has 




Moreover, we need the following lemma. 


Lemma III. 11. There exists a* G 


such that if a G then 


^{Z' 21 ) < 0, 


Vw belonging to an open neighbourhood of ojq. 

Proof. First of all recall that = 2 (tV1 — ct'^ojt, then one can compute 

^ = ^ ^_ (Vl-aJ-l)" 1 \ 

47r^/T^ \ Vl - 2aVl^^ ^ \/\F2a^- cj2 ) 

i (T^\/l + 2(T\/1 — ^ — (Vi — <7^ — l)^\/l — 2(T\/1 — cr^ 

^ A-K^ (T2(2f72 - 1) ■ 
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Since k is purely imaginary with positive imaginary part and Lrp^2P^J is self-adjoint, for the first 
summand in the expression for ‘Si{Z 2 i) one gets 

JN2{q^, g^*) \ 

- 


^{qa,A) = -25R 


^Lj,^2P^JN2{q<s,q<si *) 


Hence, 

»(Zi,) = -2» j , 

By direct computations one has 


020 ( 2 :) = (Lt — (2z^t + 0)) ^-^20 = A- 


\/u}^^+2^\x\ j \ 


47r|x| 


+ C^ 


3—*\/—‘^T-|-2^Tkl ( 1 


— I 


47r|x| 


with 


A — -^[{{2a -I- Pj^/lPr — W—ojt + 2^r)(-f^2o)i + (^\/^ + V—^t + 2CT)(-f^2o)2] 

C = ^[((20- -I- 1)^/ut — V^T + 2^t)(-H2o)i — — iV^T + 2^t)(-H2o)2] 

where d = 2i(2cr-|-l)cj 7 ’-|- 2 (fT -\-l)y^ujj'y/~^T + 2^j'—2i(cj -\-l)y^uj'j''\/ooT + 2^j'~2 \/lot 2^^^/—lot + 2^j'. 

From which follows 


47r 

9a2o = -J 


{iy/—UJT + 2 ^t — V^T + 2^t)(-H2o)i 

((2cr -I- l)-ywr — Vo’T + 2 ^t — ^/—^T + 2^r)(-H2o)i 


-h 




—i{2y/uJT + x/wt + 2^t + i^/—iOT + 2^r)(-H2o)2 
{ — y/ojT + 2^t -I- iyJ—OJT + 2^t)(-H2o)2 


Hence 


^(( 9020 ) 1 ) = [^(-f^2o)i(—(o -I- 1 )\A^\/o:t + 2^tV—‘^t + 2^r + ((o -I- l)a;T + Ct)\/— wr + 2 ^t)+ 

(-H2 o)2(—(2(ct -I- l)^r -I- (2(T -|- l)a;T)\/wT + (Ct + (2cr -|- l)wT)\/wr + 2Cr)] 

^(( 9020 ) 2 ) = [*(-f^2o)i((2(cr -t- l)^W'r -I- + 2^r — (3 ct -|- 2),/uJtV‘^t + 2^tV—^t + 2^t)+ 

+(-f^2o)2((o' -|--|- ((cr -|- l)a;'r — ^t)V^t + 2^t)]- 

(47) 


Moreover, by (27) one gets 


JA^2(g^) 


^1 


-2cj|ga;Tl (g'^)l(g®)2 

|2^-i(3(gvi/)f + (gvi>)i) 2cr((T - l)|g,^j,p'^-^ 




(48) 
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which implies 


-2ia\q^^\^--^ (l - (l + 


2a\q, 


(jJT \ 


2(T—1 f ^ Vl —( 7 ^ — 1 


H20 = {S'Ky/ujTP^JN 2 {qm)Gi^^,Gujj,)L' 2 ' = 


( 49 ) 




IGttAw^/^ V a 


+ 


+ 


^ i -(■'" 2(«»))2 (v 5 =fe^ 


\/l —(T^ —1 


O' V^OIT’ +^T + y/^T 


{JN 2 {q^))l (V^t+ 5 t+x/^ 


Let us notice that (48) and (49) imply 


i{H 2 o)ii{JN 2 {q^))i = - 2 ^(JA^ 2 (g^))f+ 


+ 


■y/^ 


y/l — a^ — 1 


AtvIk y y/LOT ~^T~\~^/^T ^ \/aJx’H“^r + \/^T 

{,H2^)2i{JN2{q^))x = (JiV2(gvl>))2i(JiV2(9v,))l + 

\/^ / 1 _I Vl-o-^-l 1 \ ( T]\J.(m')P 

AniK y -^ujx — ^ ^ ^ ^ 

^(-f^ 2 o)l(«/-^ 2 (<?'I '))2 = 25 ^*(-^A' 2 (g'>Ir))i( JA 2 (giIr ))2 + 


i(JA^ 2 ((?^))l(JA^ 2((?^))2 


+ 


\/^ 


1 ^- ^'{JN 2 {q^))l 


^TTIK, Y \/ LOt — ^T~\~-\/^T ^ \/ LOt H“^T ~^yj 

{H 2 o) 2 {JN 2 {q ^,))2 = {JN 2 {q^))l + 


+ 


y/^ _ 

AiTvitZ ^ y/LO t — ^T~\~y/ 


1 _|_ Vl-a ‘^-1 1 


O' \/^T~^^T~^y/^T 


iiJN 2 {q^MJN 2 iq^)) 2 , 


then by (46) and (47) it follows 


3f?(Z'i) = -23f? (± ^oJN2{q^) \ ^ l(3fJ((g^^J,)i(JiV2(g^))i +cj((g„^j2)(JAr2(g^))2) = 

\ Av / tK 






(T 


with 


f{a) = ^(^2(1 + cr)^ + 2ay/ 1 — Y^—1 + 4(j\/l — ^ — {2+2>a)\l —1 + 4(T'\/i — <7^1 + 4o''\/l — cj^+ 
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+ 1 + 


-i + Vi^ 


a- 


a 


^(—1 — cj)\/—1 + — 16(7^ + + cr + 2a\/ 1 — a/”! + 4(T\/ 1 — ^ 


1 


— l+'s/l—cr^ 


-l + 2a 


l+\/l —2(T\/1 —^ 0'+cr\/l+2cr\/r— 

( —l+Vl —0-2)^ 


V 


+ 


\/l —2(T-\/1 —(T^\/i+2(j\/ 1 —y 


+ ^l+cr+^l + (T — 2(T\/ 1 — 0-2^ Y^l + 4(T\/ 1 — 0-2+ 


+ 1 + 


-i + VT^ 




(T 


—1 — 2(7 + (—4(7 — 4(7^) \/l — a"^ + (1 + 2cr + 2(7v/1 — a'^j \J 1 + Aay /1 — (7" 


/ 


1 - 


+ 


— l+\/l —c 


l+\/1 —2(T\/T—(T^ (T+rr^J 1+2(T\/1 —(j^ 




V 


'\/1—2cr-\/l —cr^'\/ l+2cr'\/T^(T^ y 

Notice that one has /(a) ^/>0, and m ^ —oo as a ^ l/\/2; this implies 

l.m l.m i = o-. 

(T-5-1/V2 7r|dp (T^-i/v^ 

Hence there is a neighborhood of where 3^(Z2]^) is strictly negative. A Mathematica plot of the 
function f{a) in the range given in figure 

Summing up, one can conclude that there exists a* G 

r 1 


III.IO 


VxA’ 


(7 


j_ Vs+i 
VV2’ 2 V 2 


such that ^{Z 2 ^) < 0 for cr G 

□ 


Remark III.12. The following reformulation of the equation for zi will be useful. First of all, if we 
denote Kt = K{ujt), then the ordinary differential equation for zi becomes 


zi = iiz\_ + zATrl^^il^^i + Zn, 


for some remainder Zn. 

Secondly, let us notice that zi is oscillating while y = \zi\‘^ decreases at infinity. Hence, it is easier 
to deal with the variable y, which satisfies the equation 

y = 2ifi{iKT)y^ + Yr, (50) 


where Yr is some suitable remainder. 
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FIG. 1. f{a) 


Remark III.13. From Lemma III. gl we have 

\{g + P^hi + P^ki, Ziq)\ < TZ((jj){\\g\\L°° + — 

U) ^ U) ^ w ^ 

< — uj\ ||/i||L°°_ + + WP^kiWr^^ ), 

W W w ^ 

hence 

|Yr| = \Zr\\z\ = \ZR + i{K-KT)\zi\^zi\\z\ < 

< 7^l(a;, |z| + \\f\\L-_J\z\[{\z\^ + + \z\\ut - a;|(|z|2 + \\hU^_J+ 

+ \z \ llPj-fci ||2,oo_ P \z\ ||Pj-/ii ]• 

w ^ w ^ 


IV. MAJORANTS 

In this section we exploit the so-called majorant method to prove large time asymptotic for 
the solutions of the modulation equations. Preliminarily, we need some assumptions on the initial 
conditions. 

A. Initial conditions 

Let us fix some e > 0 to be chosen later in order to obtain a uniform control in the estimates. 
Then we assume that 
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( 51 ) 


|^(0)| < ei/2 

ll/(0)||Li<ceV2, 


where c > 0 is some positive constant. 


From the definition of zi (see Proposition III.5) one has 


zi — z = TZ{ijj)\z\^ 


Then the following estimate holds 

y(0) = |zi(0)|2 < |z(0)|2 + 7^(w, |z(0)|)|z(0)|3 < e + 7^(a;, |2(0)|)e3/2. 


We also need an estimate for the initial datum of the function h{t). For this purpose recall that, 
from the definition of / and h one has the decomposition h = f + {P^ — P^)f ■ Hence, 

IIM0)llLi < ll/(0)llL^. + llT-^’T)/(0)llit + KiH|u,t-^|||/(0)||i»_,, 

for some constant c > 0. 

Thanks to the former estimates, one can prove the following lemma. 


Lemma IV. 1. Let us assume conditions (51) on the initial data. Then 

|Z(0)|2 




< 


ce 


(l+t)3/2 - (l + t)3/2’ 


for all t > 0. 


Proof. Let us denote C = /q p{T)dT (the quantity p was defined in (33)). 

From the definition of the exponential and the idempotency of the projections one gets 


giCn± ^ + Pf 


Then it follows 

giC(n+-n-) ^ + P^)(n^e"*‘^ + n+ + P^) = n+e*^ + + Pf 


The lemma follows from the fact that Lt commutes with the projectors H^, the definition (34) of the 


operator Lm and the decay of the evolution of the functions Pfau, i, j = 0, 1, 2, stated in Lemma 


III.3, namely 
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= ||e^^*P^(e*^n+ + + P^)(a2o^^(0) + aiiz(0)z(0) + ao2z{of)\\L^_^ < 


|z(0)|2 ce 


'(l + i)3/2 - (l + t)3/2' 


□ 


B. Definition of the majorants 

We are now in the position to define the majorants: 

Mo(T) = max Ilot — w| 
0<t<T 

Mi{T) = max \z(t)\ 
0<t<T 

M 2 {T) = max \\P^hi{t)\\Loo_ 


( ' v 

(52) 

V1 + et y 

e \-V2 

(53) 

1 F et J 

■ g \ -3/2 



(54) 


0<t<T 

We shall use the following vector notation 

M = {Mo, Ml, M 2 ). 

Remark IV.2. From the estimates on g, ki and the definitions of the majorants follows 

\l°° = \\g F Pq^hi F Pfk + Pfki\\ijOo < 

'll) J- H) 


< TZi{uj) ( \ujt — a;| + \z\^ + 


(l+t)3/2 


\Pfhi\\Lco ] < 

'll) -f / 


< ^^7^l(a;)(M2 + eV2M2) 


From the assumptions (51) on the initial data one obtains 


y(0) < e + P{e^l‘^M)e’l‘^ < e(l + Pi{e^l‘^M)e^l‘^), 
\\h{0)\\Li < ce^l‘^Pi{e^l‘^M)e^Mo{\ F Mf + 


(55) 
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C. The equation for y 


We aim at studying the asymptotic behavior of the solution of equation (50) for the variable 


y introduced in Remark 2.18. To do that we need the following lemma which is the analogous of 
Lemma 4.1 in 


Lemma IV.3. The remainder Yji in equation (50) satisfies the estimate 

.5/2 




(1 + etfi\fiei 


{i + \M\y 


Hence, equation (50) is of the form 


y = 2^{iKT)y^ + Yr, 


(56) 


with 

ik{iKT) < 0 , 

?/(0) < eyo, 

where yo and V > 0 are some constants. Then we can apply Proposition 5.6 in [9] and get the next 
lemma. 


Lemma IV.4. Assuming the initial eondition and the source term of equation (50) as above, the 
solution y{t) is bonded as follows for any t > 0 

y(0) / . \3/2 


vit) - 


1 + 2'^{KT)yot 


<cY: ^ 

1 T ct 


where c = c{yo, Q{Kt))- 


D. The equation for Pfhi 


As a first step let us estimate the remianders in the equation (41) for hi. This is done in the next 
two lemmas. 

Lemma IV.5. The remainders Hr and can be estimated as 

/ . \ 3/2 

< 7^(eV2M) J ((1 + Mif + + \M\f), 
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and 


WtHUlI < 7^(eV2M) j ((1 + Mi) 3 + eV2(i + \M\fy 


Proof. From the estimate (37) on Hr one has 


WPt^rWl]^ < 77.2(0;, |z| + ||/||l“_J[| 2;|3 + (|z| + \u)T - o;|)(|zp + + 

iJJ .til 'll! 


+ ||P^/ii||loo_ ) + (|2;|2 + ||£,oo^ + ||P^/ii||i:oo_ )2] < 

lit J- in m 




1 et 


3/2 


Mf+ 


+ 


\ 1/2 \ / / ^ 3/2 

Yht) "■ + rTl«“") (+ Jiiw‘ + (iT^ 


M 2 + 




+ 


(l + t)3/2 Vl + ei 


3/2 


M 2 


< 


< 77(e3/2M) 


1 H” 6t 


3/2 


((l + Mi)3 + e3/2(l + |M|)^). 


The bound for Hf^ follows in the same way from the estimate (38). 


□ 


In the next lemma we get a estimate the evolution under the linear operator Lt of the remainder 
P^Hr. 

Lemma IV. 6. For any t, s >0 the following estimate holds 


\e^^^P^HR{s)\\L^_^ (1 + t)3/2 < 77(e^/2M) 


1 + es 


3/2 


(Mi3 + e3/2(l + |M|)3). 


Proof. From the analytic expression (42) of PIr and the estimates of the evolution of the functions 


02 O) Oil) and ao 2 stated in Lemma III.3 one has 




< 772 ( 0 ;, |z| + ||/||loo_ )| 2 ;|[| 2 ;||o;r -o;| + (| 2 ;| + + ||/ii)^] < 

in in in 
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/ , \l/2 [/ \3/2 

<n{e^/^M){— - Ml —- M 0 M 1 + 

Vl + es/ Vl + es/ 




\ 3/2 

< {Mf+ e^/^{l + \M\f). 


From the two previous lemmas we can get the following result. 


Lemma IV.7. Let us consider the equation for Pfhi 


dPfhi 

dt 


with initial condition and source terms satisfying 


, u ] — QL^{Pfhi,v) + {PfHji, v)i2 + {PfH'fi, qvGx)i2, 


such that 


11^1(0) 11 ri < e^^^ho, 


Hr = Hr + H'fi, 


WPfHRUr^ < 


\PfH'f,\\m < H 2 


^ / e 


1 + et 


_ _ / f \3/2 

\\e^^^PfHR{s)\\L^ ^ (1 + tfG < Hs {Mf + e^/\l + \M\f). 

\ 1 + es/ 

for some positive constant ho, Hi, H 2 and Ho. Then its solution is bounded as follows 


\Pfhi\\L^_^<cl,^] iho + Hi + H2 + H3), 


where c = c{ujt) > 0. 
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Proof. By the Duhamel representation one has 


{P^hi,v)L2 = ^ ^ 




for all V € V. 


Then from the dispersive estimate in Theorem VI.3 and the estimates on the remainders proved 
above in the duality paring defined by the inner product one has 


\Pfhi\\L^ = sup 


{Pf,hi,v)L 2 


< 


O^v&Li \\V\\vnLl 


(l + ()3/2ll'‘l(")lla+/ (l + /_,)3/2 


(||^T-f^R('S)llLJ„ + \\PTHRis)\\Ll)ds+ 




< c{u}t) 


3/2 


ho + 


1 P et J Jq (1 P t — 

The lemma follows from the fact that 

/■* 1 




10 


(1 +1 — s )^/^ \1 + 


es 


3/2 


ds < c 


1 T ct 


3/2 


for some constant c > 0. 


□ 


E. Uniform bounds for the majorants 


To prove that the majorants are uniformly bounded, the following lemma will be useful. 


Lemma IV.8. For any T > 0 the majorants Mo, Mi, and M 2 satisfy the following inequalities 

Mo{T) < 7^(eV2M)[(l + Mif + e(l + \M\)\ 

(Mi(r))2 < 7^(eV2M)[l + eV2(i + \M\ f], 

M2{T) < n{e^/^M)[{l P Mif P + |M|)^]. 
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Proof. It follows form Lemma IV.4 and IV.7 as Lemma 4.6 in [18], but we give the proof for sake of 
completeness. 

Step 1. Let us begin noting that 

+ ll/lli°°_i < '^2(w, |t| + ||/||L°°_J(|tP + \\PTki\\L°°_^ + \\PThi\\L°°_ ) < 

w ^ w w w 

s (irtij575 + ^ 


< 77(e^/2M)^j-^(l + Mf + 


Then by the definition of Mq and the bound on Icot — a;|: 

e 

1 T et y \ 


Mo(r) < max 
0<t<T 


(r^) (1 + Miirfp 


+e^/^M2{T)fdT + ( —^ ) (1 + Mf + 


\ ct 


< 


< TZ{e^/^M)[{l + Mi)^ + e(l + \M\f]. 

Step 2. Since y = |tip, we can exploit the inequality proved in Lemma 
77(e^/^M)(l + |M|)® and y(0) < eyo, one gets 

\ 3/2 


IV.4 


the fact that Y = 


y < 


l + etj ^ 


From which follows 


< y + TZ{uj)\z\^ < 




\ ct 
Step 3. Recall that 




3/2 


Mf 


< 77(e^/2M)[l + + |M|)' 




Hi = TZ{e^/‘^M){{l + Mif + + \M\)^), 
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H2 = n{e^/^M){{l + Mif + + \M\f), 


H 3 = n{e^/^M){Mf + + \M\f). 

Hence from Lemma llV.71 follows 

/ , \ 3/2 

WP^hiU^ , < nieP^M) —— ((1 + Mi)3 + eP\l + \M\f), 

which implies the inequality for M 2 . □ 

We are now in the position to prove the uniform boundedness of the majorants. 

Proposition IV.9. If e > 0 is sufficiently small, there exist a positive constant M independent of 
T and e such that 


\M{T)\ < M, 


for all T > 0. 

Proof. From the previous lemma follows 

\Mf < 7^(e^/2M)[(l + Ml)® + e^P{l + |M|)®] < + e^/^F{M)), 

where in the last inequality we have replaced the estimate for Mf, and F[M) is a suitable polynomial 
function. 

Furthermore, M(0) is small and M{T) is a continuous function. Hence it follows that |M| is bounded 
independent of e <C 1. □ 

The last proposition gives a summary of the behavior of the functions w(t), z{t), Pfhi(t), and f{t). 

Corollary IV. 10. There exists a finite limit ojoo for the function uj{t) as t ^ +00. Moreover the 
following holds for all t > 0 

l^oo - w(t)| < Mj^, 

\zit)\ < M , 

||/(t)||L», <M^. 
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V. LARGE TIME BEHAVIOR OF THE SOLUTION AND SCATTERING ASYMPTOTICS 


A. Large time behavior of the solution of equation Q 


The results of the previous section lead to the following theorem. 

Theorem V.l. Let u{t) he a solution of equation Q with initial datum uq ^ V H of the form 

uo{x) = + zo^{x) + zq'^*{x) + fo{x), 

where 9 q G M, wq > 0, zq G C with 

|2(0)| < ll/olli. < cfW, 

for some e, c > 0. Then, provided e is small enough, there exist uj{t), ^{t), z{t) G C^([0,+oo)) 


solutions of the modulation equations (21)-(23), and two constants lOoo, M > 0 such that Uoo = 


lim oj{t) and for all t >0 

t—¥ ~ 1“00 


i{t,x) = + zm{t,x) + zm*{t,x) + f{t,x)) , 


where 


|woo - w(t)| < 


1/2 




\ ct 


Proof. Let us recall that the decomposition of the function / as 


f = g + hi + k + ki 


depends on the quantity uj{T). On the other hand Corollary IV.10 claims that the function oj{t) 
converges to some oioo > 0 as t —)■ +oo. 

As a consequence, one can reformulate the decomposition by choosing T = +oo. Moreover, all the 
estimates obtained before for finite T can be extended to T = +oo without modification. Hence the 
theorem. □ 

The next goal is to construct precise asymptotic expressions for w(t), 7 (t), and z{t). For later 


convenience let us define (recall that ^ depends explicitly on ui, see (12); and similarly for K, see 


(44) and subsequent, and 7 ) 


f,00 - 
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Too — t(^oo) 


Kqq - K (cc^o 


Lemma V.2. Under the assumption of Theorem V.l the functions uj{t), 7 (t), and z{t) have the 
following asymptotic behavior as t ^ +oo; 


uj{t) = Woo + i ; cos(2^ooi + oi log(l + ek^t) + 02 ) + 0{t 

1 + ekoot 1 + ekoot 


Q2 


lit) = Too + log(l + ek^t) + 0{t ^), 


zit) 


e*/o 

/ , ^ 1 — 2(5 

(1 + ekoot) 2 


+ 0(t ^), 


where 


^+00 ^ 1—25 

^oo = 2;i(0)+ / + eA:oos)^Zi(s)ds, 

Jo 

ekoo = 2^iKoo)yo, d = , and qi, q 2 , ai, 02 , hi are constants. 

Proof. We will prove the asymptotics for z{t) only; the formulas for a;(t) and 7 (t) can be deduced 
as in Sections 6.1 and 6.2 of [9]. 

In order to do that let us recall that the equation for ziit) can be written as 


zi = iizi + iKoo\zi\^zi + Zr, 


moreover Remark III. 13 and the inequalities satisfied by the majorants in Lemma IV.8 justify the 
following estimates on Zr 

\Zr\ < l^l + ||/||l°°_^)[(|2|^ + ||/||l°°_j)^ + \z\\u:t — + ||h||L°°_ )+ 

7 (t J- tn ^ lit 


+ |2:| liR'p/ui ||x,oo_ “1“ 1 2 I ] T 

70 -I 111 -1 


< TZ(e^/^M)^ 


(1 + et)3/2-\/et 
as t —?• + 00 . On the other hand, Lemma |IV.4| implies 


(1 + M’) = 0(t-"), 


hit) = 


y(o) 


1 + 2Q(Koo)yiO)t 


+ 0(t-3/2), as t —)• + 00 . 
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Let us note that \zi\ satisfies the same bound of j^l, namely 

_ . ^ .1/2 
\zi\ < M 


\ €t 


then the equation for zi{t) can be rewritten in the formulas 

= t^zi + ^Koo-—7^z:pr^-z—7^zi + Zi, 

1 + 23'(iLoo)y(0)t 

where Zi = as t —)• +oo. 

Since y(0) = eyo, one has eKooVo = |e/soo(l — iS) and the equation for zi{t) becomes 


ii = (i{ - 


zi + Zi. 


Hence, one gets 

e^Jo 


zi{t) = —- OLf^i(0)+ f e + ekoos)^ ds) = Zoo 

(1 + ckoot) 2 V dO / 


e*/o 


(1 + ekoot) 

where Zoo is as in the statement of the lemma and 

/>+oo 


(1 + ekoot) 2* 


ZR 


ZR = 


pifs 


1-i S 

l + ekoos\ 2 


1 + ekoot 


Zi{s)ds. 


The bound on Zi implies zr = 0{t i). Therefore z{t) has the asymptotic behavior as t 
stated in the lemma because 


z{t) = zi{t) + 0{t = z, 


+or')- 


(1 +ekoot) 2 ' 


B. Scattering asymptotics 

Let us make the following ansatz 

u{t, x) = s{t, x) + C{t, x) + f{t, x), 

where 

s{t,x) = 

is the modulated soliton and 

C(t,x) = e*®W[(z(t) + z(t))Ti(x) +f(2:(t) - z(t))T2(x)] 


+00 


□ 
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is the fluctuating component. Recall that the functions 'hi and ^2 satisfy 








Therefore from equation 0 one gets 

0^’^) 2 ^ - \qs\‘^''qs - aiq{z+z)^i - <^ 2 q{z-z)^ 2 )^+ 

+ (7('S + C) - j“('5 + C) - + f^'2) + iz-i^z){'^i - i'I>2)],v)L2, 

for all V € V, where Qo is the quadratic form of the free Laplacian. Hence, as in [T], the solution 
f{t) can be formally expressed as 

f{t,x) = Ut* fo{x)+ i f Ut-T{x)qf{T)dT - i f Ut-r * G{T)dT, 

Jo Jo 


where we have denoted 


G{t) = + C(t)) - iuj{t) — isit) + C(t))+ 


-ie*®W[(i(t) - + i'^ 2 {t)) + {z{t) - iCz{t)){^i{t) - i^ 2 it))] 


:|2 

and Ut{x) = is the propagator of the free Laplacian 

In order to prove the asymptotic stability result we need the two following lemmas. 


m 


Lemma V.3. If the assumptions of Theorem V.l hold true, then 

rt 


'0 


p+oo ^+00 

Ut-T{x)qf{T)dT = Ut* / U-r{x)qf{T)dT - / Ut-T{x)qf{T)dT = Ut * 4>o + ro, 

Jo Jt 


where G and ro = 0{t as t ^ +00 in 


Proof. The strategy is similar to the one exploited in the case without eigenvalues (see the proof of 
Theorem 7.1 in my- since foix) = (^^,^( 372'('0 (H^)’ some function 4>o '■ K''~ —^ C, one gets 


2 

U 


(47r)2 


f 


p+oo 

V?. I 

0 


f2o{j 


r‘^dr = 


( 271)2 


r+°° ^ 

/ \(l2o{y)\‘^3/ydy. 

Jo 
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Hence (j)Q G if and only if (/>o G ^fydy). On the other hand, one can make the change 

of variables n = ^ in the integral that defines the function (j)Q and get 

^o(y) = ^ f-l Vudu, 


then cj)Q = ^qj Moreover, by Corollary 


IV.IO 


u \u 


one has 


2 

1,2 


^+oo 2 

Jo 


Qf 


y/udu < C 


n+OO 


uj Jq {u + e)‘ 

for some constant C > 0, hence the Plancherel identity implies 

'2/Tn,+ 


^du < C, 


4>o e L^{R'^,^/ydy). 


In the same way, for any t > 0 the following holds 


Fo|Il2 = 


(2vr)2 


1 / 1 
qf(t + 


u 


u 


<C 

L'^{R+,yfudu) V 1 + 


for some constant C > 0 independent of t. Which concludes the proof. 

The analogous result for the integral function Ut-r * G{T)dT requires different tools. 


□ 


Lemma V.4. Assume that the assumptions of Theorem V.l hold true, then 

ct /*+oo /*H-oo 

Ut-r*G{T)dT = Ut* I 

to 

where (fi G and n = 0{t~^G'j as t ^ +oo in 


^+oo r-\-oo 

Ut-T * G{T)dT = Ut * / U-r * G{T)dT — Ut * / U-r * G{T)dT = Ut * 4>l + ri, 

Jo Jt 


Proof. We exploit the idea used in |18j to prove Lemma 5.5. 
Step 1: restriction to the leading terms. 


From the expansions (29), (30) and (31) for d;(t), 'y{t) and i(t) —if^z{t) follow that the function G{t) 
is made by a quadratic part consisting in the terms multiplied or e®®^*)|2;ooP, with 




\/l + ekoot 

which are of order t~^ and a remainder of order t~^G^ The convergence and the decay of the 
remainder is trivial from the unitarity of Ut- Furthermore, from the analytic definition of G it 
follows that it is a complex linear combination of functions of the form 


— e ^ ^ , (X — uJqq, T z/oo,cJoo ^oo* 
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Hence it suffices to prove the lemma for the functions n(t)(5(x), where n(t) is one between 
and e®®(*^|zooP- 

Step 2: decomposition of Ut*Q. 

Let us note that we can rewrite the convolution product as follows 


kr 


Ut*Q = 


(47rit)3/2 J^3 


f ■ i^iV) C*” 4t 

/ e-*Vg(y)dy + 

Jm 


(47rit)3/2 J^3 


i 


; (^,y) , 


e * 2 t (e* 4t - l)Q{y)dy = 


•N_ -NI 

^ \ 4t ^ f X 


(57) 


where gt(y) = (e* 4* _ i)Q(y). 

Since |e®^ — 1| < 0 and the function G(?/) is exponentially decaying as \y\ —>■ +oo, the norm of the 


second term of (57) can be estimated in the following way for any f > 1 

1 


(2t)3/2 




2t 




for some constant C > 0. Hence, recalling that n(r) < (1 + ekooT) we obtain 

^+oo 


/*+oo 

/ Il{T)Ur*QtdT ^ 

Jo 


and 


r*+oo 


U{T)Ur * Qtdr = 0{t 


-1-1 


as t —)■ +00 in L‘^{R^). 


Step 3: Analysis of the first term in (57) in a particular case. 


Let us first show how to treat the terms with the phase 0(t) replaced by cOoot. 
Note that 


Q{x) = 


1 


I I 19’ 

a + |a:|^ 

Hence, in the case of the summands with |tooP it suffices to prove the integrability of the function 


l{x) = f 
Jo 


OO I |2 

•{ 2 ; \ 


(1 + ekooT){\x\‘^ + 4ar2) 


dr = 


r 

A{x) 

Jo 


•/ 1^1 ^ 


4a 


T\/t \ 


(1 + e/coot) efeoo (kp + 4ar2) y 


dT+ 
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+ 


4a 


r 

A{x) / 
Jo 






(|xp + 4ar2) 


and the decay of 


Itix) = e 


OO I |2 

•/ \^\ \ 
^(a;ooT-■^)_ 


T 


(1 + ekooT){\x\‘^ + 4ar2) 


dr = 


Mx) Jt 


“ l£l!) / Vt 


•t hC \ 

oH^^ooT-n^) 


4a Ta/t 


V(l + e/coor) e/coo (Ixp + 4ar2)) 


dr+ 


+ 





.Ml 

4t 


(|xp + 4aT2) 


dr 


/i,t(x) + /2,t(a;), 


where A(x) = 4 ^ 4 t||xp • 

For the function l 2 {x) one has 


\h{x)\ < 


4a 

e^/c^ 


r 

A{x) 

Jo 


p r2,^3 

(|x|2 + 4aT2) 


With the same estimate it is trivial to prove 


h,t{x) = 0{t 


as t —)■ + 00 , in L^( 

In order to treat Ii note that 

4a Ta/t 


Since 


+ 


(l + ekoot) e/coo (|x|2 + 4aT2) ekoo^/T{l + ekoox) e/cooV^(|xP + 4ar2) 

= as t —)■ + 00 , is integrable on (0, + 00 ) and A{x) G L^(]R^) one has to 


ekaa^{l+ek 

00 


prove 


\x\‘^A{x) 


^+00 

Jo ' 


■/ 1^1 ^ 


-v/r(|xp + dar^) 


dr = 


r +00 

A{x) 

Jo 


I 1^1^ ^ 1 


4cjooT -C?7" — 


^+00 

4a4.(x) 

Jo 




r3/2 


/o \/^ ^0 

From formulas 3.871.3 and 3.871.4 in [16] one has 


(|xp + dar^) 


dr G L 


2/in>3 


p+00 |^|2 1 pi-nA 

A{x) / —^dr = -G 

do V'T 
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It remains to handle with the second integral in the former sum which can be done integrating by 
parts in the following way 


A{x) / J . —^dr 

Jo 


(|xP + dar^) 


= AA{x) 


>0 


dr 


r7/2 


(|rcp + 4ar2)(|xp + dwooT^) 


dr 


< 


r+oo 

<CA{x) 

Jo 


r5/2 




0 (|x|2+ 4min{a,a;oo}r2)2^^ ^^ 

Then we are done. 

In order to estimate the decay of Ii^t it suffices to study the decay of 

e/coo Jt ^/T{\x\'^ + Aar'^) 


r+oo . ,2 1 f+oo , ,2 3/2 

= A{x) / _ 4aA{x) / ^ 

' Vt Jt (|rc|2 + 4ar2) 


dr, 


Jt V' Jt 

which can be done integrating by parts as before. Let us do that for the second term (the computation 
for the hrst one are analogous and simpler): 

Aix) / gh^oor-^) r 

Jt (|x|2 + 4ar2) 


= AA{x) 


r+oo 


i(a;oor-hL)_(^ 


dr 


r7/2 


(|x|2 + 4ar2)(|x|2 + AojoqT^) 


dr 


< 


< CA{x) 


t-l/2 + 


r*+C>0 


r5/2 


't (|x|2 + 4min{a,a;oo}r2)2 


dr 


< 




The case of the summands with is analogous, while the case of Zoo^ is more difficult because 
|xp + 4(a;oo - 2 ^ 00 )+^ = 0 for 


T = t* = 


2^ 

00 ^00 
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Let g : M"*" —>■ M"*" be a continuous function with the properties: 


0 < g{t*) < t* W > 0, and A{x)g{t*) G 


It follows that g{t*) = 0{t*) = 0(|x|) as |x| —>■ +oo. Hence, one can represent (0,+oo) = (0, 
g{t*)]U {t* — g{t*),t* + g{t*)]U {t* + g{t*), +oo). Integrating by parts once more one has 


A{x) 




g*((t*^oo 2 ^oo)t 4 ^ )_ 


^3/2 


\x\^ + 4ar2 


dr 


< 


< CA{x) {t* - g{t*))-^/^ + 




^5/2 


(|x|2 +4aT2)||a:|2 + 4(woo - 2^oo)r^ 


-dT+ 




+ 


^9/2 


(|x|2 + 4ar2)2||xp + 4(a;oo - 2^oo)r^ 


■dT+ 




+ 


^9/2 


(|xp + 4ar2)||x|2 + 4{uJoc - 2^oo)t2|2 


dr < 


< CA{x){{t* - + (f - g{t*)fh e 


where the last inequality follows from formula 3.194.1 in [16]. In the same way (exploiting formula 
3.194.2 instead of 3.194.1 in [IHI)) one has 


L 


>o I |2 f3/2 

gi((i^oo—2^oo)'r )_ 

+git*) |x|2 + 4ar2 


dr 


< 


< CA{x){{t* + g{t*))-^/^ + {t* + g{t*))-^/^ + {t* - g{t*))-^/^) G L‘^{R^). 


Finally, 


L 


t*+g{t*) 


2 ^ 00 )^— 141 ^) _ 


^ 3/2 


t*-git*) 


IxP + 4ar2 


dr 


< 


< CA(x) [ ^ -^dr < G l2(m3). 


Summing up, the integrability of the integral function 

/•H-OO 


/o 


U{T)Ur * Qdr 
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is achieved. It is left to study the decay of 


A{x) J e 


+00 , ,2 ^3/2 

,i((a;oo-25oo)T-LL) * 


\x\‘^ + 4ar^ 

First of all, let us note that integrating by parts one obtains 


dr. 


*((<^00 2 ^ oo)t )_ ~ _ 

(O,t*-g(t*)]n[i,+oo) |xp + 4ar2 


A{X) I g^VV-oo--soo;.-^ 


^3/2 


< 


< CMx) 

< CA(x) 


d 


^7/2 


dr (|x|2 + 4ar2)(|xP + 4(a;oo - 2 ^oo)t2) 


dr < 


\x\‘^ + ^ J t ||x|2+ 4(a;oo - 2 Coo)t 






^9/2 


+ 


zdr . 


(|x|2 + 4aT2)||x|2 + A{uJoo - 2Coo)t2|2 
The three integrals in the last inequality can be estimated in the following way: 


(i) ir"' ’ s c/f” S 


t 

f\ rt*-g{t*) 
(ii) Jt 




_ r**-g{t*) 




|a;P+4(2^oo— lJoo)t 2 /i (|a;|+2v^2^oo-t<JooT)||x|-2v^2^oo—(.^oor| 


dr 


< , dr < cr 1/2; 

•-^0 ||a?| 2\/2^oo ^oo'^l 

iiiii < 0-1/2 

^ ^ Jt (|a;P+4ar^)||xP+4(a;(x,—2^oo)r^P — JO 

< C't- 1/^(1 + In ||x| - 2^6000 - 2 .^ 00 (t* - fl'(T'))|). 


dr 


Hence, since y4(x)ln||x| — 2\/a;oo — — git*))\ £ one can conclude 

A{x 


' (0,t*—5'(t*)]n[t,-|-oo) 

as f —)• + 00 , in 

Let us now observe that 


i((<^cx> 2^cxd)t ^= 0(t 1/^) 


|x|2 + 4ar2 


H(x) / e 

J {t* —g{t*),-^oo)r\[t^-\-oo) 


"^((^oo 2,^00)'?' LI. )_ 


^3/2 


|x|2 + 4ar2 


dr 


< 


< CA{x) 


r+oo 


d 


f/2 


dr (|x|2 + 4ar2)(|x|2 + 4(woo - 


dr < 


dT+ 
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< B{x)A{x) ^ — 2^'^ — CB{x)A{x)t G L^(M^), 

where i? : —>■ M"*" is a continuous bounded function. 

Finally, 


A{x) 


^i((u}ao 2Joo)t LI. ) 


t3/2 


'{t*-g{t*),{t*+g{t*)]n[t,+oo) 


IxP + dar^ 


dr 


< 


rt*+g{t*) 1 

< CA{x) / < C'A(x) 5 (r)t -^/2 g 

di V"^ 


Summing up, thanks to the unitarity of Ut, we proved 

/ + 00 

U{T)Ur * Qdr = 

as t —)> + 00 , in L^(]R^). 

Step 4: conclusion of the proof. 

The conclusions of the previous step hold true if the phase Uoot is replaced by &{t). In fact, the 
estimates which involve the integral of the absolute value are totally unaffected by change of phase, 
then it is only left to adjust the argument involving integration by parts. This can be done integrating 
by parts exactly as before, which leaves a factor g*(©{*)-‘^oob in the integrand. Then, the boundary 
terms can be treated in the same way because |e*(0(*)-^°°*)| = i. Finally, the extra contribution to 
the integrand can be estimated as it is done for the summand arising from differentiation of 
since |0(t) — Wool < t t > 0, where C is a positive constant. □ 

Summing up, we have proved the following asymptotic stability result. 


Theorem V.5. Let o ^ , for a certain a* € and u{t) € C(R~^,V) be a solution 

of equation 0 with 

u(0) = no = +^)^i + i{zo - ^)^2] + h^Vn 


for some loq > 0, 70 , zq G M and fo G n 

no is close to a solitary wave, i.e. 


Furthermore, assume that the initial datum 


l^^ol 


< and Wfoh},. < 


where c, e > 0. 
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Then, if e is sufficiently small, the solution u{t) can be asymptotically decomposed as follows 


u{t) = +Ut*foo + roo(i), as t ^ +oo, 

where Uoo, ek^o > 0 , 6i G M and foo, ?’oo(i) £ with 

\\roc{t)\\L 2 = 0{t~^^^) as t +00, 
in 


Remark V.6. Numerical evidences (see Lemma III.11} suggest cr* = — 0,96. 


VI. APPENDICES 


In the following appendices we collect auxiliary material on the model studied. In Appendices A 
and B we recall results from [3] regarding resolvent, spectrum and dispersive behavior of linearization 
operator L (see ©) to help the independent reading of the present paper. In the subsequent 
appendices C,D,E we state and prove further details regarding spectral properties of the linearized 
operators, in particular the structure of eigenfunctions associated to the discrete spectrum, the 
structure of the generalized eignevectors, and the proof of Lemma II.8. 


A. Resolvent and spectrum of linearization 

We denote 




47r|x| 


ca ^ 0, A G 


—XGx'2* 


47r 

1 

> 

co . 

M 

to 

_1 

1- 

> 

to 

* 

— A^A2* 

^ IT(A2)' 

-*Si A 2 


(58) 


with the prescription ± iX > 0. 

Furthermore, we make use of the notation {g, h) := g{x)h{x) dx. 

The resolvent of linearized operator L is described in the following 

Theorem VI.1. The resolvent R{X) = {L — XI)~^ of the operator L is given by 
R{X) = 

where 

IT(A^) = 327r2aia2 — 4i7r(ai + 02 ) (^V—oj + iX + V—w — iX^ — 2y/—oj + iXy/—oj — iX, 


( 59 ) 
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and formula (59) holds for aZZ A G C \ {A G C : W{\^) = 0, or 3f?(A) = 0 and |9‘(A)| > ui}. 


Furthermore, the symbol * in (59) denotes the convolution and 


= 7r*T {Gu,-ix{x) - G,^+ix{x)) , r;)^2(x) = \ {G^-i\{x) + G^+iA(a:)) • 

ziX 2 

Finally, the entries of the second matrix are finite rank operators whose action on f £ L^(M^) reads 
Aif = [iA(47rQ;2 - W-oj + i\){Qx 2 , f) - (47rai - ix/-w + iA){Tx 2 , f)\G^+ix+ (60) 

+[iA(47ra2 - iV-w - i\){Qx 2 ,f) + {^nai - ixf-u - i\){Tx2, f)]Gui-ix, 

A 2 / = [iA(47rQ!i -W-u + iX){Qx2,f) - (47ra2 - iV-w + iA)(rA2, /)]G^+a+ 

+[iA(47rai - W-uj - iX){Qx 2 ,f) + (47ra2 - W-u - iX){Tx 2 , f)]G^-ix, 

El/ = -[zA(47ra2 - W-u: + iX){Qx 2 , f) - (47rai - ix/-oj + iX){Tx 2 , f)]G,^j^ix+ 

+[iA(47ra2 - iV-w - iX){Qx 2 ,f) + - ix/-uj - iX){Tx 2 , f)]G^-ix, 

E 2 / = -[iA(47rai - W-uj + iX){Gx' 2 , f) - (47ra2 - ix/-oj + i\){Tx 2 , f)]Gt^+ix+ 


+[iA(47rai - iV-uj - iX){Qx2,f) + (47ra2 - iV-w - iX){Tx2, f)]Gu,-ix- 


The resolvent determines also the spectrum of the operator L, as given in the following 
Proposition VI.2. The spectrum of the linearized operator L has the following structure 

(a) (Tess(T) = {A G C : 5?(A) = 0 an(i|3‘(A)| > w} 

(b) If a £ (0, llx/2), the only eigenvalue of L is 0 with algebraic multiplicity 2. 

(c) If (T = lly/2, L has resonances hiuj at the border of the essential spectrum and the eigenvalue 
0 with algebraic multiplicity 2. 

(d) If a £ { 1 /^/ 2 ,1), L has two simple eigenvalues ±i^ == Iii2ax/l — 0 ‘^uj and the eigenvalue 0 
with algebraic multiplicity 2. 

(e) If a = 1, the only eigenvalue of L is 0 with algebraic multiplicity 4. 

(f) If a £ (1, + 00 ), L has two simple eigenvalues ±2aV <7^ — loj and the eigenvalue 0 with algebraic 
multiplicity 2. 
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B. Dispersive estimates 


We recall that the propagator is the inverse Laplace transform of the resolvent. So the 
dispersive behaviour associated to the linearized dynamics projected on the continuous spectrum is 
controlled by the following result, proved in [3], Theorem 4.8. 


Theorem VI.3. Let ^ 


There exists a constant C > 0 such that 


1 

27ri 



{R{X + 0) - R{X - 0)){x)e-^^f{y) dXdy 



for any f G where 


C+ = {A G C : 51?(A) = 0 and ^(A) > w}, C_ = {A G C : 3?(A) = 0 and 9(A) < —w} . 


C. Eigenfunctions associated to and generalized eigenfunctions 

1. The eigenfunctions associated to 

Here we describe the eigenspaces associated to the simple purely imaginary eigenvalues = 
±i2cjVl — 

Let us start with the eigenvalue The following proposition holds true. 

Proposition VI.4. The eigenspace associated to if is spanned by 



Proof. In order to prove the proposition we need to solve the equation 

LT = 

in D{L). For x 7 ^ 0, the previous equation is equivalent to the system 

f (-A + = 0 

\ ^2 = |(-A+ w)4'i 

from which follows that 'Ll must belong to and solve the equation 

(-A + UJ - C)(-A + UJ + O'Li = 0. 
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Hence, the solutions in 


are of the form 
^i(x) = A 




4tt\x\ 


4it\x\ 

^-^u)+e,\x\ 


^2(x) = iA ^ - iB 

for any A, B G C. 

It is left to look for A, B G C such that 'I'i G D{Li) for i = 1, 2, i.e. 

-^A - ^B = -{2a + 1)^{A + B) 

-i^A + i^B = -^{iA - iB) 

Exploiting the fact that ^ = 2ay/l — a'^uj one can show that the two equations of the previous system 
are linearly dependent and 

\/l - 0-2 + 1 , 


B = -- 


a 


The thesis follows by setting A = 1. 

Let us note that in the previous proof we have chosen the constant in such a way that ^i(x) G 
and 'L 2 (x) G for any x G \ {0}. This fact will be used to prove the next proposition. 

Proposition VI. 5. The eigenspace associated to —i^ is spanned by 


□ 


= 


^'1 




Proof. In the previous proposition we proved that 

L2^2 = 

^ -LiTi = iC'i’2 

with Ti real and T 2 purely imaginary. 

Taking the conjugate of both equations and recalling that the operators Li, i = 1, 2 act on the real 
and imaginary parts separately, one has 

L2(—'I'2) = — 

-LiTi = -if{-^2) 


which is equivalent to 


LT* = 


because the operators Li, i = 1, 2 are linear. The proof is complete. 


□ 
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2. The generalized eigenfunctions 


Our goal is to compute the generalized eigenfunctions associated to the continuous spectrum. In 
order to do that, we treat the two branches C+ and C- of the continuous spectrum separately. 

Proposition VI. 6. The generalized eigenfunctions associated to are 


-^u+ri\x\ / 1 \ -iv^ri=U|x| / 1 \ pi^/u^\x\ / 1 


for any rj £ [u, +oo) and D £ C, with 

A= . ^ , {C + D), 

^ _ {2(j+l)uj+{a+l)^{i^/rP^—y£j:fLj)—iy/rtTTf2 ^ 

— {2a+l)u)+{(T+l)y£o{i^/iT£I:+y/rj+uj)—iy/rp^-oT 
Proof. For any rj £ [oj, +oo), we need to solve the system 




where G L°°(M^) does not necessary belongs to As in the computation for the eigen¬ 

function at if X 7 ^ 0 the former equation is equivalent to the system 

f {-A + uj-f,){-A + (jj + f,){^+)i = 0 
\ (^+)2 = |(-A + a;)(^'+)i 

which leads to 


^'+(x) = A- 


,-^LU+r]\x\ 


p^LO+r]\x\ / 1 \ / 1 \ pWV-^\M I 1 


for some A, B, C, D £ C. Since we require 'h+ G we get B = 0. Moreover, the bound¬ 

ary conditions in the domain of the operators Li and L 2 must be satisfied by ('I'+)i and ('I '+)2 
respectively. Then A, C, and D solve the system 

I -^A - i^C + i^D = - + C + D) 

I ,VAMA+^C-^D = -^{-iA + iC + iD) 


which concludes the proof. 


In the same way, one can prove the analogous result about C_. 
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Proposition VI. 7. The generalized eigenfunctions associated to C- are 


-^uj-r]\x\ / 1 \ -^^/-('l^+u])\x\ 

^-(x)=A . , , I \+C- 


47r|x| 


47r|x| 


1 \ / 1 

+ D- 


—i 


47r|x| 


—i 


for any rj G (—oo, — w], where D G C and 


A = , "fl.' , (C + D), 

p-, _ (2cr+l)aj+((T+l)v^(i'^-(r?+aj)-^ 

-(2a+l)uj+(a+l)A^(i-^-(ri+oj)+y'r;+oj)-iyt'r)^-uj^ 


It is easy to see that the projection operators from 


onto X^, X^ and X^^ are given by 




I »(/■>■■•) ,I,., 

pcf = f_ pOf _ plf^ 


duj I 

K = 0(4',4'*), 


Moreover, we denote with 11^ the projections onto the branches C± of the continuous spectrum 
separately. 


D. Proof of Lemma III.8I 


In this appendix we prove Lemma II.8 whose statement is recalled for the reader’s convenience. 


Lemma VI. 8. There exists a constant C > 0 such that for each h G X^ holds 

||[PV-i(n+-n-)]h||^^<c||h||i»_^. 

Proof. From the definitions of the operators P'^ and 11^ one gets 

pv - f(n+ - n-) = n+(j - u) + n~(j + u) = 


1 

27rf 


' (R(A + 0)-R(A-0))(J-iI)dA+ / (P(A + 0) - P(A - 0))( J + f/)dA 

Jc+ Jc- 


We will estimate just the hrst integral because the second one can be handled in the same way. 


Exploiting the explicit form of the resolvent (VI.l) it follows that 
R{A){J-iI) = {iAgp*+Tp*) 


1 i 

dvr 

Ai + S 2 

i(Ai + S 2 ) 

-i 1 

^ L»(A2) 

—i{M + El) 

A 2 + El 
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= i?*(A)(J - il) + RmiX)iJ - il), 


where iZ* and Rm correspond to the convolution term of the resolvent and the multiplicative term. 
Note that 


{x - y) + Tx 2 {x - y) = 2Guj-i\{x - y) = 


-uj+i\\x-y\ 


2tt\x - y\ 

is continuous on C+. Hence, the integral on C+ of the convolution addends vanishes. 

Let us now consider the multiplicative addends in the integral on . From the explicit formulas for 


Ai and S 2 given in Proposition VI. 1 one can compute 

(Ai + S 2 )(x, y) = 87r(a2 - ai)G^_ix{y)G^+ix{x) + [87r(a2 + ai) - W-w - iX]G^_ix{y)G^_ix{x) = 


-u)+iX(\x\ + \y\) 


Denote 


Then it follows 


D±(A2) = D((A±0)2). 


[(.I?m(A + 0) — Rm{X — 0))(J — iI)]i^idX — 


Jc+ 

f + {{a + l)^/io - W-u - i\)e^V-‘^+iXi\x\+\y\) 

lc+ 7r|x||y|D+(A2) 

f a e^x/-ux-i\\x\ ((g, i\/-w - 


dA+ 


Jc+ 'K\x\\y\D-[X^) 

If we compute the change of variable k = xj—u — iX in the first integral of the last equality, and 
k = —xf—oj — iX in the second one, then one has 


[(.I^m(A + 0) — Rm{X — 0))(J — iI)]i^idX 


Ai 


'K\x\\y\ 


lc+ 


r'+oo „-\/P+2 L|j/| /■+00 p-y/kA+2u)(\x\ + \y\) 

<yV^k -e-'^l"ldA:+ / 2ik^ 


D{k) 


< G- 


=-\/AA|i/| 


min|p,e-^l^l| < G 


D(k) 


,-VXiX\y\ -VXk}\x\ 


dk 


< 


\y\\x\ [|x|’ j \y\ |a:| 

where the first inequality is obtained integrating by parts both integrals. 

The integral of the other three elements of the matrix operator (i?m(A + 0) — Rm{X — 0))( J — il) can 
be estimated in the same way and this implies the statement of the lemma. □ 
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